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In order to provide experimentalists with a useful tool for experiments on the fc-boson nonlinear 
coherent states, we investigate some critical features of these states, as regards their realization as 
, stationary states of a trapped laser-driven ion. After illustrating the fc-boson case, we analyze the 

■ single-boson nonlinear coherent states and prove that their coherence and minimum-uncertainty 
properties are violated whenever the Lamb-Dicke parameter is one of the roots of the Laguerre 

£\J , polynomials which enter such states. We discuss the possibility of circumventing this problem 

, , . through the relaxing of the rotating wave approximation. 

PACS numbers: 03.65.-w, 42.50.Dv, 03.65.Fd 

o : 

. I- INTRODUCTION 

Coherent states have become a tool of the utmost importance in quantum optics and many efforts have been devoted 
fS |' to investigate them, also with respect to their prospective applications in the growing field of quantum technologies 
(QT) such as quantum information (coding, transmission and elaboration of information by exploiting properties of 
£3 ' quantum states), quantum computing, quantum imaging and quantum metrology. 

In particular, their original definition has been generalized, so that nowadays several quantum optical extensions 
of the Glauber states have been constructed and analyzed from both the theoretical [H]-[II| and the experimental 
P om t of view. Such states can find application also in quantum calculus, as shown, e.g., in 
In the present paper the focus is on the fc-boson and single-boson nonlinear coherent states (NLCS) which can 
be realized experimentally by means of trapped ions interacting with driving electromagnetic fields, as reported in 
references Q and [lj|, respectively. Due to their mathematical structure these states are affected, for specific values 
of the Lamb-Dicke parameter, by singularity problems which do not appear to have been tackled systematically so 
ly-^ far. The paper is organized as follows: in section 2 the Hamiltonian of a ion trap driven by two lasers is illustrated, 
while the construction of the corresponding fc-boson NLCS is reported in section 3 resorting to standard tools of 
quantum optics, i.e., the interaction representation and the rotating wave approximation (RWA). Sections 4 and 5 
are devoted to the special case of the single-boson NLCS. The relaxing of the RWA in order to remove the singularity 
' problems is discussed in section 6. 

o : 

• • . II. THE PHYSICAL SYSTEM 

> • 

■ The properties of the ion-trap systems and of their Hamiltonians in various physical regimes have been widely 
investigated [2(|-[29j]. As for the construction of coherent states for ion-trap systems, in [7] it has been proved that 
single-boson NLCS can be realized by exploiting the Hamiltonian H = Hq + H mt of a two-level ion trapped in an 
external harmonic oscillator potential and interacting with two laser fields. 

The free Hamiltonian Hq describes the motion of both the electronic (or internal) and vibrational (or external) 
degrees of freedom of the ion, while the interaction Hamiltonian H iat corresponds to the interaction of the ion with 
the electromagnetic fields of the driving lasers. The external states of the ion are the ground state \g) and the excited 
state |e), their transition frequency is u) = u) e — LJ g , v is the frequency of the trap harmonic potential. With a the 
annihilation operator for the vibrational motion of the trapped ion, h = a^a, a z the Pauli matrix, which acts as a 
third component pseudospin operator, 

Hq = hvh + — Juva z . (1) 

In Eq. (p} the zero energy level of the ion has been set at halfway between \g) and |e), so that in outer product 
notation a z — |e)(e| - |<7)(<?|- The interaction Hamiltonian reads 

H int = A (e e^ox-uot) + Ei e i(k lX - Ul t) | s+ + H c _ 
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In Eq. ([2]) A is the dipole coupling matrix element, Eq, E\, and loq, w i ar e the amplitudes and the frequencies of 
the lasers fields, k , k\ their wavectors; x = r]{a^ + a)/fc^ is the position operator of the ion center of mass, where 
kh ~ ~ k\ is the value of both lasers wavevectors, r\ = kLy/h/{2Mv) is the Lamb-Dicke parameter, with M the 
ion mass. With a x and a y Pauli matrices, 5+ = [a x + icx y )/2 = \e){g\ and S- = S\_ are the two-level electronic flip 
operators. For the lasers tunings one sets luq — ^, which corresponds to the resonant condition of the first laser, and 
uji = u) — kv, which means that the second laser is tuned on fc-th lower (or red-sideband) vibrational level, k > 1. 
From Eq. © one then obtains 



S 4 



H.c. 



(3) 



III. THE fc-BOSON NLCS 



The underlying idea for the construction of nonlinear bosonic coherent states from the Hamiltonian ([3]) can be 
summarized as follows. 

1. Through the unitary transformation 

U(t) = cxp (— iHot/h) = exp (—ivht) exp {—iuia z t/2) , 

one writes the Hamiltonian © in the interaction representation Ti m t = (t)Hi nt U(t); utilizing the RWA "Hint 
is put in the form [191 ] 



(4) 



^oc(f< > S + + H.c. 

where is a non-Hermitian operator acting in the Fock space; 

2. the condition Ff®\£k) = [f|, [III determines the states upon expanding them in terms of number states; 
besides, (£fc|?4nt |£k) — by construction; 

3. the states are recognized to be the eigenstates of the annihilation operator 

A k = f(h)a k , 



(5) 



with the appropriate identification of the nonlinear function fin), i.e., they are coherent states in the Barut- 
Girardello sense [30| . 

Let us detail the above steps. In the interaction representation one obtains 

Eq ^ ( ir l) m+e r ivt( m -i) a j m a e 



Hint = XEie^i 12 



E x E m\l\ B 

£,m=0 



E 

e,m=o 



V-V) %vt(m-l+k) t m t 

f v 'a' a 



S+ + H.c. 



(6) 



The link between the physical system of the laser-driven ion trap and its mathematical structure is more clearly 
grasped rewriting Eq. ^ in the form 



E, 



+ Q k 



S+ + H.c. 



(7) 



where O and Qk are the operators 



^ (ft + 1) • • • (h + r) (h + 1) • • • (n + r) n w ; 



iurt rH/Jl 



o = 



(8) 
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{h +1)-.. (ri + s + fc) ^ 



+ e**4°V) + £ at s " fc 6 ^ ^ . (9) 

^-^ (n + 1) • • • (n + s — k) 

s=k+l \ ' I \ ' ) 

In Eqs. ©, ©, L^\r] 2 ), with a integer > 0, is the 'La guerre polynomial operator'. By this we means an appropriate 
integral representation [3l| . Its action in the Fock space results into the generalized Laguerre polynomial of degree n 
and parameter a 32] 

n + a \ x 



= E(-)'("„t;) 

£=0 V / 



fl ' (10) 



Applying the RWA to the interaction Hamiltonian Q implies that all terms oscillating with frequency ^ or larger 
than v are disregarded in Eqs. (jHJ and ([9]) so that in this approximation O and reduce to 



6(0) =4°) , Qi 0) = {iri) k J[{n + 3)- l Lf{r?)a k 



i=i 



This leads to the form ((4|) 



where 



^ = \E 1 e-^/ 2 Pl 0) S+ + H.c. , (11) 



P (0) = ^ d (0 )+ Q(0)^ (12) 

The states which are the stationary states of the master equation of the system, are found from the condition 
Fjf ] \( k ) = 0. Their physical meaning is twofold: they are both noiseless codes in the sense specified in (33[, as 
(£fc|?4nt = 0, and the coherent states of the system, defined as the eigenstates of the operator ([5]) 



where 



while the nonlinear operator /(n) is 



4t|&> (13) 



k j (fc) / 2\ 

/w^n^+.r 1 ^. d5) 

oo 

( 0) IA.\ = n with IA.\ = A k )\ n ) 



One finds that the following self-splitting into k states results from the condition \£k) — with |^) = c 

n=0 

fe-l 

l&) = E • ( 16 ) 



«=0 

The states |£fc)^, 



n=0 
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are fixed by the corresponding sets of coefficients 

„(*) _ ^3 



-nk+e 



y/(nk + £)l 



n-l 



n w + mk ) 



(18) 



(k) 

c) being the normalization constants, chosen to be real to retrieve the Glauber states when k = 1 and / = 1 



£ 

n=0 



(nk + £)\ 



n ^ + mfc ) 



(19) 



In Eqs. (TIE)) , (|f 9p the scalar functions f(£ + mk) are obtained from the action of the operator (|15p in the Fock space, 
n-i 

while Y\ = 1 for n = 0. 

m=0 

It is readily verified that the states (fTT)) are eigenstates of the operator (0, Ak\£,k)e = £,k\£,k)e, with the fc-degenerate 
eigenvalue (fT4")) ; since j(^fc|Cfe)^ = <5j'.£j they form an orthonormal basis of a fc-dimensional Hilbert space. On the other 
hand, due to the splitting (fl~6)) . \£ k ) is eigenstate of Ak with the nondegenerate eigenvalue k£k- 



IV. THE SINGLE-BOSON NLCS 



In the nondegenerate case k = 1 one retrieves the construction reported in Q ■ The notation adopted so far can be 
simplified dropping the index k throughout. 

A. The states |£) as eigenstates of A — f(n)a 

Eqs. (fl~3|) and (fT4|) show that |£) is by construction eigenstate of the operator §5§ with k = 1, A = f(n)a, the 
eigenvalue being £ = iEo/(r)Ei). The relevant commutation relations are [n, .4] = -A, [n, A f ] = A* and 



[A, A f ] = (n + l)/ 2 (n) - n/ 2 (n - 1) = C(n) . 



(20) 



In terms of number states 



10 =Mj2c n \n), 



(21) 



where the normalization constant TV, which is the new label of in Eq. (|19p . has been factorized so that from Eq. 



< " " 1 n_1 L (0) fr) 2 l 



(22) 



3=0 



with co = 1, f(j) = Lf\rj 2 )/[(j + 1) LyV)] and 



TV = 



£m 



,n=0 



(23) 



With m > 1, the following recurrence formula for the coefficients (|22|) holds 

-l 



c 



(n + m) ! 



n /( n+ ^) 



(24) 
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B. The states |£) as minimum-uncertainty states 



For the dimensionless operators X = (A> +A)/V2, P = i(A f - A)/y/2, in view of Eq. ([2"0|). one has [X,P] = iC(n), 
while the variances 

(AX) 2 = X - [(AA) 2 + (AA^) 2 + 2{A*A) - 2\(A)\ 2 + (C(n))] , (25) 

(AP) 2 = -\ [(AA) 2 + (AAt) 2 - 2(A*A) + 2\(A)\ 2 - (C(h))] , (26) 
calculated in the coherent states minimize the Heisenberg uncertainty relation (AX) 2 (AP) 2 > \([X, P])| 2 /4 

(AX) 2 = (AP)f = \ mC(nM)\ ■ (27) 

OO OO 

InEq. © (Z\C(n)\0 =M 2 Y. \cn\ 2 C(n) = ^p n C(n), where C(n) = (n + l)f 2 {n) - nf 2 (n - 1) and Pn = \c n \ 2 M 2 

n=0 n=0 

with < p n < 1 Vfi. Eq. (f27| confirms that the states (|21l) are intelligent states al ong with the Glauber states. On 
the other hand they are not coherent states according to the Perelomov definition [3J| since, in view of (I2TJ1) . A and 
At generate an infinite-dimensional algebra. 

V. THE SINGULARITIES 



Eq. 



shows that for n > 1 c„ has n(n — l)/2 zeros and as many poles due to the roots of the n Laguerre 



polynomials Lj (r/ ) and Lj (if), respectively, whose effects on the coherent states (f^Tj) have to be investigated. 

It is worth acknowledging that the above problems are referred to briefly in [23| where it is stated that the values 
of rj 2 must be such as to prevent them. Furthermore, in [24j the effects of an ill-defined function f(j) are investigated 
in relation to the representation of the Fock states in terms of single-boson NLCS. 

A. A single zero of c„ 

Let us investigate the effect produced by a single zero, rj 2 = fj 2 , of a given c n , say cp. This means that L { ° ] _ t (f) = 
and, equivalently, f((3 — 1) — > oo. Since cp — 0, from Eq. (|24|) it follows that cg +m = 0, V m > 1. Due to the above 
conditions the sums in Eqs. (|21[) and (|23[) can be truncated so that the states and their normalization constant read 



(28) 



n=0 



■p-i 

Ei 

,n=0 



(29) 



The vanishing of the coefficient cp, violates the coherence of the states (|28|) as the latter do not satisfy the eigenvalue 
equation of the annihilation operator 

Notice that the truncactions (|28|) and (|29|) do not remove the singularity of the states which influences the variances 
of the observables X, P. In fact, addressing now the minimum-uncertainty property in the states \£)p, from 

(AA) 2 = e [uj^^-Ml^f-Mj^Y 
and from Eqs. (J25j), (p6]l. one gets (AX) 2 = A + E, (AP) 2 = -A + 3, where 

a = l(t 2 +e 2 ) [pp-m-pj,-i(p)-pp-2(P)] 

p-i 



lep^m [i - pp-itM + \ Y,pnW)c( 



n=0 
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FIG. 1: The product (AX)l (AP)j for /3 = 4, |£| = 0.1 



while p„(/3) = | Cn | 2 AAj , < p ra (/3) < 1 Vn, /3 and * denotes complex conjugation. 

As an example, in Fig. Q]the product (AX)1 (AP)2 = S 2 — A 2 is plotted in the range < i] 2 < 5 assuming that 

the singularity problem corresponds to j3 = 4, i.e., r? 2 is one of the roots of L^\t] 2 ). Of course, in Fig. Q] r/ 2 varies 
with continuity so that the product of the variances contain the roots in the considered range of all the Laguerre 
polynomials involved in the /3 = 4 case, i.e., L^(r) 2 ) and Ly(r] 2 ), j — 1,2,3. One emphasizes that a root of the 
polynomials L^\rj 2 ) (the single zero case) corresponds to a singularity of the product (AX) 2 (AP) 2 , while when a 

root of the polynomials Lj(rf) occurs (the single pole case of Sect. 5.2) the product vanishes. 

Notice also that in the example considered the minimum- uncertainty product 1/4 is still maintained over a wide 
range of values of ry 2 . This feature strongly depends on the parameter /3, in the sense that the lower the value of 
P, the wider the ranges where, due to the minimum-uncertainty property, the truncated states (|28[) can still be 
considered, approximately, intelligent states and conveniently utilized for QT purposes. 

B. A single pole of c n 



Here rf is a single pole of a given c M , i.e., c p — > oo. This implies that il_i(?? 2 ) = and correspondingly f(p — 1) 
: 0. From the recurrence relation (I24[) it follows that c^ m — > oo, V m > 1, as rapidly as c^. The above conditions 



imply that the contribution of the first /x terms to the sums ([211) and (|23|) can be disregarded with respect to the 
singular terms, leading to the following expressions of the states and their normalization constant 



\0m = K X] c « 



(30) 



E 



(31) 



Extracting the singular coefficient c M from both Eqs. (|3U|) and (|3"Tj) and utilizing the recurrence property ([2"4")) . one 
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finds 



10, 



1 pn m 

[so*)]-* E 4i n /~ x wi 



e =f i 



where is the series 



so*) = E 

j=0 



i=0 v¥+f) 1 <=/i 



I*' n / W 



(32) 



The question whether or not the normalization constant S(/j,) converges, and hence the states ([3U[) exist (in the sense 
that they are normalizable) , can be tackled, for example, through the ratio test. According to this criterion the 

oo 

convergence of the series ([32p . written as S(fi) = E] ' •> implies the evaluation of the ratio 

3=0 



2+1 







(m + j + i)/ 2 (m + .?) 



r(l) 
L »+3 



(v 2 ) 



when j — > oo. Resorting to the asymptotic formulas for the Laguerre polynomials of large degree reported, e.g., in 
35], one finds that 



lim 



S 3 + l 
b 3 



j-s-oo 



tan [2<J(p + j)?f 



(33) 



The result ([33]) proves to be wildly oscillating as j increases, which means that the very definition of the states is 
questionable when a pole of a coefficient c M is encountered. This pathology becomes manifest in the example of Fig. 
Q] where the poles correspond to the vanishing of (AX)| (AP)|. 

The above analysis shows that there are values of rj 2 which alter significantly the structure of the single-boson 
NLCS, resulting either in the truncated states (|28| or in the states (|30| whose definition is not certain. These 
singularities are inherent in the very construction of such states which relies on the application of the RWA to 
the two-level ion model. One is therefore led to investigate whether relaxing the RWA can remove the above difficulties. 



VI. RELAXING THE RWA 



A detailed analysis of the possible corrections of the RWA is out of the aims of the present paper. Here we discuss 
two ways through which a 'first-order' correction of the RWA can be implemented, where first-order means that only 
the terms of Eqs. ([8]), ([9]) associated to the exponentials exp(±ii4) are retained. Specifically, one adds the following 
time-dependent operator ?{S to the interaction Hamiltonian ([lip 



where 



= XEl e-^/ 2 



Ei 



Sa 



H.c. 



= (it}) 



n + 1 n + 1 



(34) 



(fc+i) 



2^ e -ivt 



,k+l 



(h + 1) • • • (n + k + 1) 



+ 5 M e- fc '4 V) 



(irj) 



k-i r (fc-i) 



k-l 



(h + 1) • • • (n + k - 1) 



(35) 
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The first approach develops along the general scheme of section 3, i.e., one assumes that, in this approximation, 
the coherent states \z), 



\ z ) = ^2 d n\n) 



n=0 



are determined from the condition P\z) = 0, t = F 1 (0) + PP with J F 1 (0) from Eq. Jl2|) and = (Eq/E^ (1) +q[ 1 \ 
where and are given in Eqs. and (l3"5j) . The coefficients d n can be calculated from the relations 

■^-7= d 2 + irje- wt w(t)di + w(t)d = , (36) 

x /(n + l)(n + 2) V« + l 

+«;(t)L<°V)d» + ^^i^K-i = , n > 1 , (37) 

where = Eq/E\ +exp(ivt), if <io an d di are known. To hx the latter coefficients one could resort, e.g., to suitable 
orthonormality conditions. 

The second approach entails the utilization of the RWA according to the procedure proposed in Q: one defines 



J? 00 / . \2H t (1) / 2\ 

= ^0 \^M^ a tV + i V $ iff j a = A (0) - SoM, so that = -£ /£i|0- The condition J» = 

1 fci V* 1 ) n + 1 

with .F = G^ + PP, is then imposed upon expanding \z) in terms of the coherent states |£) 

\z) = I > 



on account of the completeness relation / d/i(£)|£}(£| = I, d/i(£) = d 2 £/7r, with d 2 £ = d(Re£)d(Im£), being the 

00 

measure in the complex £-plane. Since P[ \£) = A n (£;t)\ri), where the coefficients A n (£;t) are readily calculated 



n=0 

in view of Eqs. fl2T}, ([21]), one has 



j» = / d/i(o<€k> 



/CO 

n=0 
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In Eq. (39, the expansion of |n> in terms of the coherent states \?), \n) = /<MO«»IO, with <£» = <(£'), 
has been utilized. Left-multiplying F|z) = by J d/i(£")(£"| and resorting to Eq. (|3"8"|) . one finds 



z = i^+jd^')m,e-,t), (39) 

where /C(£, £) denotes the kernel 

lC(Z,e-,t) = X)A n (€;tK(0- 

n=0 

The time-dependent shift ([39| of the eigenvalue —Eq/Ei = ir]£ of the operator G^ implies that also the values of 
the Lamb-Dicke parameter, possibly giving rise to a singularity, are shifted from the pathological value to a close 
one, most likely not singular. 
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VII. CONCLUSIONS 



To summarize, we have analyzed the properties of the NLCS, which correspond to the experimental scheme of 
laser-driven trapped ions. Upon verifying that they can be defined as both Barut-Girardello states (eigenstates of 
annihilation operator) and minimum-uncertainty states, we have carried out a mathematical investigation whose 
main conclusions, in view of physical applications, are: 

• For single-boson NLCS, if the squared Lamb-Dicke parameter happens to be a root of a given Laguerre polyno- 
mial i^_ 1 (?7 2 ) both the coherence and the minimum- uncertainty properties do not remain valid; however, the 
latter feature is still maintained over a wide range of typical values of r\ if /3 is not exceedingly high (i.e., one 
does not consider very excited vibrational levels); if the squared Lamb-Dicke parameter happens to be a root of 

a given Laguerre polynomial L^\(r) 2 ), the very existence of the single-boson NLCS is compromised, due to the 
oscillating behaviour of the defining series; 

• the general construction which leads to the fc-boson NLCS is naturally affected by the same kind of pathologies; 

• we have proposed to relax the RWA conditions to get around the above difficulties and have discussed, for the 
single-boson NLCS, first-order corrections of the known theory; 

• from an experimental point of view our results suggest that, due to the finite accuracy in fixing the parameter 77 
when preparing single-boson and fc-boson NLCS for actual experiments, one should be cautiously far from the 
singularities analyzed in this paper in order to maintain the RWA approach. 
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